Flexible generation of correlated photon pairs in different frequency ranges 
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The feasibility to generate correlated photon pairs at variable frequencies is investigated. For 
this purpose, we consider the interaction of an off-resonant laser field with a two-level system pos- 
sessing broken inversion symmetry. We show that the system generates non-classical photon pairs 
exhibiting strong intensity-intensity correlations. The intensity of the applied laser tunes the degree 
of correlation while the detuning controls the frequency of one of the photons which can be in the 
THz-domain. Furthermore, we observe the violation of a Cauchy-Schwarz inequality characterizing 
these photons. 
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I. INTRODUCTION 

The nature of light has always intrigued mankind and 
its study has nowadays culminated in the field of i 
turn optics investigating matter-field interaction 
With the first measurement of an intensity-intensity cor- 
relation function by Hanbury Brown and Twiss |3| and 
the theoretical basis for the characterization of light by 
Glauber [3, Hi, scientists have efficient tools in their hands 
to probe light fields for quantum signatures @. In the 
last few decades, the interest in non-classical light has 
grown significantly with the advent of quantum computa- 
tion and information science Q • Entangled photon pairs 
turn out to be indispensable for many quantum proto- 
cols @ and quantum algorithms [t| . Currently, there are 
a series of experimental techniques available to produce 
entangled photons such as parametric down conversion 
(10l - [l2l ] , four- wave mixin g I13l 4l5| , electromagnetically in- 
duced transparency [H, or cavity QED [18|, EH • Fur- 
thermore, an atomic memory for correlated photon states 
has been realized experimentally, playing an essential role 
for quantum communication over long distances [2(il - l22j |. 
Recently, a heralded entanglement source of great prac- 
tical importance has been demonstrated 23|, 24 1. In ad- 
dition, theoretical considerations have predicted the gen- 
eration of a correlated photon pair in the x-ray regime 
from strongly driven atomic ensembles |25|. Very re- 
cently, a communication network for quantum informa- 
tion processing has been proposed (2(| , which consists of 
numerous different nodes and channels. Since such differ- 
ent nodes may have different characteristic frequencies, 
there is great interest in investigating non-classical pairs 
of photons of different frequencies [27| . As an important 
milestone in this direction, entangled photons of different 
but close frequencies limited to the microwave or optical 
ran ges have been generated and detected experimentally 

Based on this background, we investigate here a two- 
level system with broken inversion symmetry which is 



driven by an off-resonant laser field. By means of adjust- 
ing the laser frequency wi, one can spontaneously gener- 
ate a photon at an approximate frequency u>l — ljq and 
a subsequent photon with transition frequency ujq. With 
the parameters of e.g. gamma globulin macromolecules, 
those frequencies can be in the THz- and optical regimes, 
respectively, see Fig. [T] We find that this photon-pair of 
different frequencies is both of non-classical character and 
entangled because it violates a Cauchy-Schwarz inequal- 
ity. The advantage of our scheme lies in the fact that the 
frequency of the longer wavelength photon can be manip- 
ulated by an appropriate selected detuning. This is quite 
useful in driving a quantum network composed of differ- 
ent nodes of various frequencies including quantum wells 
or dots of THz transition frequencies. Furthermore, the 
high flexibility distinguishes our model from a cascade 
three-level system or other down conversion processes. 



II. THE MODEL 

In particular, we consider a two-level system (see Fig. [T]) 
with the transition frequency ujq described by the or- 
thonormal ground state |1) and excited state |2) with 



Bloch Siegert shift 
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FIG. 1: (color online) The emission of the non-classical pho- 
ton pair. The non-resonant laser excites the two-level system 
with broken inversion symmetry and induces the emission of 
a THz-photon and the subsequent spontaneously emitted op- 
tical photon. 
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broken inversion symmetry, meaning that the diagonal 
parts of the dipole operator satisfy the following con- 
dition: | 7^ | g>22 1 : where we define pij — e{i\r\j) 
for {i,j} E {1,2}. The system is driven by a classi- 
cal off-resonant laser field given by a linearly polarized 
monochromatic plane-wave field in the dipole approx- 
imation E = Eq cos(wit) with laser frequency u>e and 
amplitude Eq . The sample is surrounded by a quantized 
environment that accounts for the processes of sponta- 
neous emission [30| . 

The Hamiltonian H describing the system takes into ac- 
count the energy of the environment He and of the two- 
level system Ht, the interaction-energy between the laser 
and the two-level system Hn and the interaction-energy 
between the environment and the two- level system H12'. 
H = H E + H T + H n + H12, or 

H = fkj k a k a k + fiLUaSz + hil(S + + S~) cos(tJit) 

k 

+HGS Z cos{uj L t) +iJ2(9k • d){a{ - a k )(S+ + S~), (1) 



The time-dependent part can be regarded as a perturba- 
tion to the time-independent part and we can thus ap- 
ply the second-order perturbation theory [H, [33| , since 
G < we, fi < ojl and (g k ■ d) < ul- 



H, 



pert = I dtH". 



(5) 



Our final Hamiltonian Hf = H' + H pert acquires the 
shape 

Hf = H(ujk - WL)a k ak + fi(^o - cjl + 

k L 

+ ^(S+ +S-) + iJ2(9k-d)(a{S- -a k S+) 



_3G_ 

8iuiL 

n 

2iujL 



^(gk-d^alS+e 1 ^ -a k S- 

k 

^2(gk ■ d)(a k - a k )S z , 



-iu) L t\ 



(6) 



where we define the Rabi frequency f2 = p\2 ■ Eq/%, 
and G = fp n — P22) ■ Eq/H leads to broken inver- 
sion symmetry [31|. Here, we have introduced the usual 
atomic operators S + = |2) (1|, S~ = |1) (2| and S z — 
(|2) (2| — |1) (l|)/2. a k and a k denote the creation and 
annihilation electromagnetic field operators of the fc-th 
mode of frequenc y uj k . The coupling constant g k is de- 
fined as g k — y/2nhaj k /Ve\, where e\ is the photon 
polarization vector, A E {1,2}, and V is the electro- 
magnetic field quantization volume. The electromag- 
netic atom-field interaction is given in the usual dipole- 
approximation. We stress the fact that we do not work 
in the rotating wave approximation, but rather choose a 
perturbative approach to account for non-linear effects. 
For this purpose, we first perform a unitary transforma- 
tion on H with Hq = J2 k hwLa k a k + HojeSz, 



H = ei Hat (H - H )e--k Hot , 



(2) 



which may be separated (H = H' + H") into a time- 
independent part 

H' = %{u) k - Lo L )a\a k + h(u - lo l )S z 



(S+ + S-)+iJ2(9k- d)(alS- - a k S+), (3) 



and a time-dependent part containing fast oscillating 
terms 

H" = — S z (e iu!Lt + e-^*) 

+ ^s+e^ + S-e- 2 ^) 
+ i^9k ■ d){alS+e 2 ^^ - a k S~ 



-2iuJLt \ 



(4) 



where we keep the slowlyest oscillating time-dependent 
terms only. We notice that the time-dependent terms 
are proportional to G and are thus important for the 
description of a system with broken inversion symme- 
try. The ratios G/uje and CI/u>e are small for optical 
frequencies lol such that higher orders are negligible in 
the Hamiltonian. Our perturbative approach also reveals 
an effect of strong driving fields - the Bloch-Siegert shift 
TiVP I (4u>l) [34| of the upper state of the two-level system, 
see Fig. [1] Finally, the two-level approximation applies 
because CI/ujl -C 1 and \ujq — wl\/ojl -C 1. 
In what follows, we shall derive the master equation em- 
ploying the Hamiltonian in Eq. © and the Heisenberg 
picture. We assume that the matter-field interaction 
is weak in the sense that an emitted photon does not 
react back on the atom and use the well-known Born- 
Markov approximation. Thus, the time-evolution of an 
arbitrary atomic operator Q(t) is governed by the Heisen- 
berg equation: ^ (Q(t)) = j- {[Hj,Q]) . By inserting the 
final Hamiltonian, we obtain 



dt 

E 



(Q(t)) 

(9k ■ d) 



[H ,Q]) 



{(a] t [S-,Q]) + ([Q,S + ]a h ) 



— ((4 [S+, Q}) + (IQ, S~]a k ) e-^) 
^{(al[S,,Q]) + ([Q,S,]a h ))}, 



(J) 



where H = h{u) - u L + CI 2 /(4u L ))S z + m(S+ + S~)/2. 
To further simplify the analytical formalism, we have to 
express the annihilation and creation operators as a func- 
tion of atomic operators in the Born-Markov approxima- 
tion. First, we insert a k (t) in the Heisenberg equation 
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and obtain the general solution for the linear inhomoge- 
neous differential equation of first order. Then we con- 
sider the leading order in the coupling and neglect the 
Lamb-shift, so that the creation operator acquires the 
shape 

4(*) = 4(0K Afct + ^(g k ■ d)S z (t)S(A k ) 



7r — - — S^(t)d(uj k - w - l^ - )) 



(8) 



where A k = u) k —LJL. We notice that for the annihilation 
operator a k , we only have to take the H.c. of the above 
formula. If we further define the different decay rates of 
the system 



1R = 7T > j ^ ( H Wfe _ W ° 

k 



7L = 

k 



K 2 

{9k • df 



n 2 



), 



h 2 



S(uj k - u) L ), 



(9a) 
(9b) 



(gfc ■ df n 2 

— d(u k - uj L + loq + J, (9c) 



and insert Eq. @ in Eq. ([7]), we may write down our 
final master equation 

j t m))-\{[Ho,Q\) 

= - 1R {(S+[S- ,Q\) + <[Q,S+]S-» 

- ^7l((S z [S-,Q]) + ([Q,S + )S z )) 

Of 

- (—) 2 it((S-[S+,Q}) + {[Q,S-]S+)) 

8uj l 

- ^-Jr({S + [S z ,Q}) + ([Q,S z }S-)) 

- (^-) 2 7l((S z [S z ,Q]) + ([Q,S Z ]S Z )), (10) 

which may be interpreted as follows: the first term ac- 
counts for the spontaneous emission at resonance ujq + 
ri 2 /(4wi), taking into account the Bloch-Siegert shift. 
The second term describes the spontaneous emission at 
the laser frequency ojl preceded by an excitation. The 
third term corresponds to the emission at frequency 
u!i — ujq — r2 2 /(4wi) preceded by an excitation of the 
two-level system. With the used parameters, later on, it 
has THz-frequency while the main resonance is optical. 
The fourth term accounts for a spontaneous emission at 
resonance preceded by an excitation (off-resonant as al- 
ways). The last term contributes to the dephasing of the 
system. We are interested in correlations between the 
processes of the first and third summands that are illus- 
trated in Fig.[T] For this purpose, we need to define these 
correlations and their time-dependent behaviors. 
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FIG. 2: The steady-state inversion operator as a function 
of the Rabi frequency with transition frequency uo = 5.0 x 



If) 1 



laser frequency ujl = ojq + 10 s , detuning A = 



10 13 s 1 and decay rate 70 = 3 x 10 6 s 1 with respect to wo- 



In order to probe the quantum nature of our generated 
photons, we calculate its intensity-intensity correlation 
function defined as [l|, Q 



(a) , . {E\-\t)E)->(t + T)EY>{t + r)Er>{t)) 



?(+), 



(+) 



(*)><£<-> (t)25f'(t)) 



?(+), 



(11) 



We know from the definition of the quantized electric field 
[l[ that E(~) oc at and that £|W oc a k . In our case, we 
also know from Eq. ([5]) that for THz-emission a\, oc S~ 
and that for optical emission a\. oc S + . Therefore, the 
probability for detecting an optical photon after a THz- 
photon as a function of atomic operators is given by 



,(2W _ (S-(t)S+(t)S-(t)S+(t)) 



Si2 J (0) = 



{S-(t)S+(t))(S+(t)S-(t)} 



(12) 



and the probability for detecting an optical photon fol- 
lowed by a THz-photon reads 



(2W = (S+(t)S-(t)S+(t)S-(t)} 
921 { ' (S+(t)S-(t)){S-(t)S+{t)}' 



(13) 



As a concrete system, we consider gamma globu- 
lin macromolecules [351 ] with the following param- 
eters |w 2 - Ui\ = 4.8 x 10 15 s~\ \p2i\ = ID and 
\p22 — pn\ — 100 D. We notice that the transition fre- 
quency is optical and we do observe the necessary bro- 
ken inversion symmetry. We choose the laser detuning 
such that the long wavelength photon is in the THz do- 
main. Alternative systems are quantum dots, which are 
0-dimensional quantum systems ha ving an electron con- 
fined in all three space dimensions [36|. Gallium nitride 
devices for example show broken inversion symmetry and 
have typical values of \p22 — Pii\ — 10 D, \pi2\ — 10 D 
and w = 4.92 x 10 15 s" 1 
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FIG. 3: The steady-state second-order intensity-intensity 
correlation function describing the probability of (a) THz- 
emission followed by an optical emission and (b) optical emis- 
sion followed by THz-emission as a function of the Rabi fre- 
quency fi. Otherwise, we use the same parameters as in Fig. [2] 



III. RESULTS 

At first, we display the population inversion (S z (t)) as a 
function of the Rabi frequency f2 in Fig. [2] We observe 
that for low Rabi frequencies f2, the population remains 
in the ground state. At a frequency of about 10 12 s~ 1 , 
we notice an increase of the population and at 10 13 s _1 , 
we see that there is a non-vanishing probability to find 
the system in the excited state. Now, we turn to the 
plot in Fig. [3Ja) of the second-order correlation function 

(2) 

$12 (0) as a function of the Rabi frequency il describing 
the probability of the emission of a THz-photon and the 
subsequent emission of an optical photon. We observe a 
strong correlation which decreases with rising Rabi fre- 
quency To induce the emission of a THz-photon, the 
system has to be excited from the ground state to the 



upper state, where it may spontaneously emit an opti- 
cal photon. Thus, at low Rabi frequencies, the emission 
of an optical photon is almost always preceded by the 
emission of a THz-photon. This explains the high de- 
gree of correlation of the photon pair. As ft increases, 
there is a non-vanishing probability to find the system in 
the excited state and an optical emission that is not pre- 
ceded by a THz-photon is possible. This means that the 
correlation decreases. Finally, we discuss the intensity- 

( 2) 

intensity correlation function g 21 '(0) in Fig.^b) describ- 
ing the probability of detecting a THz-photon right after 
an optical photon. It turns out that this probability is 
very low as expected. It slowly rises with increasing Rabi 
frequency Q. 

In this context, we also investigate the violation of the 
Cauchy-Schwarz inequality 

ff[fto)^(0) > fog^O)] 2 . (14) 

(2) (2.) 

The correlations g lx (0) and g\ 2 (0) vanish trivially and 
in Fig. IS^a, b), we notice nonvanishing cross-correlations 
violating Eq. (|14l) . Thus, we are dealing with a non- 
classical pair of correlated and entangled photons. 



IV. SUMMARY 

In summary, we have investigated the interaction of a 
two-level system with broken inversion symmetry and 
an off-resonant laser field. Using the parameters of 
e.g. gamma globulin macromolecules or certain quan- 
tum dots, we have found the possibility to generate a 
long wavelength photon in the THz-regime followed by 
a photon in the optical frequency range. Furthermore, 
we have observed a high degree of correlation between 
these photons and even a violation of a Cauchy-Schwarz 
inequality. This proves the non-classical character and 
entanglement of the photon pair. In the emerging field 
of quantum information science, non-classical correlated 
or even entangled photon pairs of different frequencies are 
of great interest, finding applications in the realization of 
a quantum network. 



[1] M. O. Scully, M. S. Zubairy, Quantum Optics, (Cam- 
bridge University Press, 1997). 

[2] D. F. Walls, and G. J. Milburn, Quantum Optics, 
(Springer - Verlag 1995). 

[3] R. Hanbury-Brown, and R. Q. Twiss, Nature (London) 
177, 27 (1956). 

[4] R. J. Glauber, Phys. Rev. 130, 2529 (1963). 

[5] R. J. Glauber, Phys. Rev. 131, 2766 (1963). 

[6] R. Loudon, Rep. Prog. Phys. 43, 913 (1980). 

[7] M. A. Nielsen, I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 



2000). 

[8] R. Horodecki, P. Horodecki, M. Horodecki and K. 

Horodecki, Rev. Mod. Phys. 81, 865 (2009). 
[9] R. Cleve, A. Ekert, C. Macchiavello, and M. Mosca, Proc. 

R. Soc. Lond. A 454, 1969, 339 (1998). 
[10] D. N. Klyshko, Photon and Nonlinear Optics, (Gordon 

and Breach Science Publishers, New York, 1989). 
[11] Han Zhang et al., Nature Photonics 5, 628 (2011). 
[12] G. Yu. Kryuchkyan, L. A. Manukyan, Phys. Rev. A 69, 

013813 (2004); H. H. Adamyan, G. Yu. Kryuchkyan, 

Phys. Rev. A 69, 053814 (2004). 



■5 



[13] S. Du, J. Wen, M. H. Rubin, and G. Y. Yin, Phys. Rev. 

Lett. 798, 053601 (2007). 
[14] M. Macovei, G.-x. Li, Phys. Rev. A 76, 023818 (2007). 
[15] M. Kiffner, M. S. Zubairy, J. Evers, and C. H. Keitel, 

Phys. Rev. A 75, 033816 (2007). 
[16] V. Balic, D. A. Braje, P. Kolchin, G. Y. Yin, and S. E. 

Harris, Phys. Rev. Lett. 94, 183601 (2005). 
[17] A. Fountoulakis, A. F. Terzis, and E. Paspalakis, Phys. 

Rev. A 73, 033811 (2006). 
[18] Y. L. Lim, Almut Beige, Phys. Rev. A 71, 062311 (2005); 

Almut Beige, W. J. Munro, P. L. Knight, Phys. Rev. A 

62, 052102 (2000). 
[19] Jie-Hui Huang, Shi-Yao Zhu, Phys. Rev. A 74, 032329 

(2006). 

[20] A. Kuzmich, W. P. Bowen, A. D. Boozer, A. Boca, C. 

W. Chou, L.-M. Duan, H. J. Kimble, Nature (London) 

423, 731 (2003). 
[21] K. Edamatsu, G. Oohata, R. Shimizu, T. Itoh, Nature 

(London) 431, 167 (2004). 
[22] C. H. van der Wal, M. D. Eisaman, A. Andre, R. L. 

Walsworth, D. F. Phillips, A. S. Zibrov, M. D. Lukin, 

Science 301, 196 (2003). 
[23] C. Wagenknecht, Che-M. Li, A. Reingruber, X.-H. Bao, 

A. Goebel, Yu-Ao Chen, Q. Zhang, K. Chen, J.-W. Pan, 

Nature Photonics 4, 549 (2010). 
[24] S. Barz, G.Cronenberg, A. Zeilinger, Ph. Walther, Na- 
ture Photonics 4, 553 (2010). 
[25] L. Jin, J. Evers, and M. Macovei, Phys. Rev. A 84, 

043812 (2011). 
[26] H. J. Kimble, Nature 453, 1023 (2008). 



[27] X. Guo, C. Xie, Y. Li, Phys. Rev. A 84, 020301(R) 
(2011). 

[28] C. Eichler, D. Bozyigit, C. Lang, M. Baur, L. Steffen, 

J. M. Fink, S. Filipp, and A. Wallraff, Phys. Rev. Lett. 

107, 113601 (2011). 
[29] M. Pelton, Ph. Marsden, D. Ljunggren, M. Tengner, A. 

Karlsson, Optics Express 12, 3573 (2004); J. Fan, A. 

Dogariu, L. J. Wang, Optics Letters 30, 1530 (2005); O. 

Alibart, J. Fulconis, G. K. L. Wong, S. G. Murdoch, W. 

J. Wadsworth, J. G. Rarity, New J. Phys. 8, 67 (2006). 
[30] M. Kiffner, M. Macovei, J. Evers, and C. H. Keitel, in 

Progress in Optics, edited by E. Wolf (North Holland, 

Amsterdam), Vol. 55, 85 (2010). 
[31] W. Yang, X. Song, Sh. Gong, Ya Cheng, Zh. Xu, Phys. 

Rev. Lett. 99, 133602 (2007). 
[32] R. Tan, G.-X. Li, and Z. Ficek, Phys. Rev. A 78, 023833 

(2008). 

[33] D. F. V. James, Fort. Phys. 48, 9, 823 (2000). 

[34] F. Bloch, A. Siegert, Phys. Rev. 57, 522 (1940). 

[35] V. A. Kovarskii, and O. B. Prepelitsa, Optics and Spec- 
troscopy 90, 351 (2001). 

[36] C. Kittel, Introduction to Solid State Physics (Wiley In- 
ternational Edition, 2005). 

[37] O. V. Kibis, G. Ya. Slepyan, S. A. Maksimenko, and A. 
Hoffmann, Phys. Rev. Lett. 102, 023601 (2009). 

[38] D. P. Williams et al, Phys. Rev. B 72 235318 (2005). 

[39] T. Bretagnon et al, Phys. Rev. B 73 113304 (2006). 

[40] P. Rinke et al, Phys. Rev. B 77, 075202 (2008). 



